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CALCULATED INTERACTION ENERGY AND DIPOLE MOMENT 
IN COLLISIONS OF TWO HYDROGEN MOLECULES 
by R. W. Patch 
Lewis Research Center 

SUMMARY 

Ab initio calculations of the interaction energy and electric dipole moment of two 
colliding hydrogen molecules were carried out on a digital computer using the ortho- 
gonalized valence -bond theory of McWeeny. Methods of calculating accurate values of 
all molecular integrals are given in some detail for cases with equal or unequal orbital 
exponents (screening constants). Linear, planar perpendicular, quadrilateral, and non- 
planar perpendicular configurations were included with intermolecular distances from 
2.5 to 5.5 bohrs (0. 13 to 0. 29 nm) and internuclear distances in the hydrogen molecules 
of 1. 401446 and 1. 450000 bohrs (0. 0741599 and 0. 0767292 nm). A weighted average of 
the interaction energies of the four equilibrium configurations was then taken for each 
intermolecular distance. 

The interaction energies obtained appear to be valid for intermolecular distances less 
than 4. 5 bohrs (0.24 nm) but are too large at appreciably greater intermolecular distances 
because of insufficient allowance for electron correlation in the wave function. A simple 
average interaction energy (Morse potential) for practical applications is cited which 
correlates experiments for intermolecular distances greater than 4. 0 bohrs (0.21 nm) 
and correlates the calculations in this report for intermolecular distances less than 
4. 0 bohrs- (0. 21 nm). 

The electric dipole moments obtained with the full valence -bond model appear to be 
approximately correct for intermolecular distances up to 4. 5 bohrs (0. 24 nm), but for 
appreciably larger distances are too small because of limitations of the model. Deriva- 
tives of the electric dipole moment with respect to the internuclear distances were also 
calculated. 



INTRODUCTION 


In high-temperature propulsion devices such as gas -core nuclear rockets, an impor- 
tant mechanism of heat transfer is radiant energy exchange between volumes of gas and 
between the gas and the wall (refs. 1 and 2). When such devices use high pressure hydro- 
gen, it is necessary to know the strength of pressure -induced infrared absorption to 
calculate the heat transfer. This strength depends on the interaction energy and electric 
dipole moment of two colliding hydrogen molecules (H 2 - H 2 ) and has not been measured 
much above room temperature. At the temperatures that occur in these high -temperature 
propulsion devices, smaller intermolecular distances become important in calculating 
the strength. 

The object of the work herein described was to calculate interaction energies and 
electric dipole moments for small intermolecular distances. The results should also be 
useful for calculating high -temperature transport properties of gases containing hydrogen 
and for calculating the equation of state of high-pressure hydrogen for application to light- 
gas guns. 

Past work on H 2 - H 2 interaction energy may be divided into ab initio calculations, 
semiempirical calculations, and experiments. These three types of work will be dis- 
cussed in subsequent paragraphs. In none of these interaction energy calculations was 
the dipole moment calculated, but the dipole moment by itself was calculated in a semi- 
empirical calculation by another investigator. 

Ab initio calculations for H 2 - H 2 have been carried out in 10 studies. For those 
unfamiliar with this type of calculation, a review of the methods used has been given by 
Slater (ref. 3). The H 2 - H 2 calculations may be divided into two groups: those where 
approximations were used for some molecular integrals, and those where all molecular 
integrals were calculated accurately. Since the interaction energy is a sensitive function 
of the molecular integrals, this is an important distinction. Griffing and Vanderslice 
(ref. 4) calculated energy of linear H^. Griffing and Macek (ref. 5) calculated the energy 
of square H^. Ruffa and Griffing (ref. 6) computed the interaction energy of linear 
H 0 - H 9 . All three of these studies used approximations for some of the molecular 
integrals. Accurate values of the integrals were used by Taylor (ref. 7) for the energy 
of linear H 4 , Magnasco and Musso (refs. 8 to 12) for the energy of quadrilateral and 
skewed H 2 - H 2 , Magnasco, Musso, and McWeeny (ref. 13) for the energy of quadrilateral 
and skewed H 2 - H 2 , and by Schwartz and Schaad (ref. 14) for the energy of linear H^. 
None of the studies using accurate values of the integrals were done for a sufficient 
number of configurations to permit proper averaging of the interaction energies over 
nuclear configuration space, although this is necessary for practical applications. 

Semiempirical calculations of H 2 - H 2 interaction energy have been made by 
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(a) Linear configuration A with equilibrium internuclear distances 
(ref. 261 for the two molecules. 
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(b) Linear configuration B with greater-than-equilibrium internuclear 
distance for molecule 1 and equilibrium internuclear distance for 
molecule 2. 
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(c) Planar perpendicular configuration C with equilibrium inter- 
nuclear distances for the two molecules. 


O 

1.45 

O 


1. 401446 
O -O 


(d) Planar perpendicular configuration D with greater-than- 
equilibrium internuclear distance for molecule 1 and equi- 
librium internuclear distance for molecule 2. 
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(e) Planar perpendicular configuration Ewith equilibrium inter- 
nuclear distance for molecule 1 and greater-than-equilibrium 
internuclear distance for molecule 2. 


o c 

? ? 

1. 401446 

1.401446 

O < 

3 O 


(f) Quadrilateral configuration F with equilibrium internuclear 
distances for both molecules. 
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(g) Quadrilateral configuration G with greater-than-equilibrium 
internuclear distance for molecule 1 and equilibrium internu- 


clear distance for molecule 2. 
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(h) Nonplanar perpendicular configuration H with equilibrium inter- 
nuclear distances for both molecules. 
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(i) Nonplanar perpendicular configuration I with greater-than- 
equilibrium internuclear distance for molecule 1 and equilibrium 
internuclear distance for molecule 2. 


Figure 1. - Nine configurations of two hydrogen molecules used in calcu- 
lations. Intermolecular distance R. had values of 2. 5, 3.5, 4.5, and 
5.5 bohrs (1 bohr * 5.29167X10" 11 m). Internuclear distances in the 
figure are in bohrs (not to scale). 



Margenau (ref. 15) (corrected by Evett and Margenau in ref. 16), Mason and Hirsch- 
felder (ref. 17), Vanders lice and Mason (ref. 18), and Abrams, Patel, and Ellison 
(ref. 19). A semiempirical calculation of the H 2 - H 2 electric dipole moment has been 
made by Van Kranendonk and Kiss (ref. 20). 

Experiments pertinent to H 2 - H 2 average interaction energy involve hydrogen vis- 
cosity and equation of state measurements. A large number of measurements of these 
two types have been made, but most were not at high enough temperatures to give any 
indication of the H 2 - H 2 interaction energy at small inter molecular distances. However, 
Guevara and Wageman (ref. 21) measured hydrogen viscosity to 2340 K. Fisher (ref. 22) 
fit Guevara and Wageman' s data together with lower temperature viscosity and equation - 
of-state data with a Morse intermolecular potential, which should therefore be valid 
down to an intermolecular distance of 4. 0 bohrs (0. 21 nm). 

In this study a full valence -bond, ab initio calculation of H 2 - H 2 interaction energy 
and electric dipole moment was carried out with accurate values for all molecular inte- 
grals. In order to cover nuclear configuration space, linear, planar perpendicular, 
quadrilateral, and nonplanar perpendicular configurations (see fig. 1) were included with 
intermolecular distances from 2.5 to 5.5 bohr (0. 13 to 0. 29 nm). In four of the config- 
urations (fig. 1(a), (c), (f), and (h)) both H 2 molecules had equilibrium internuclear 
distances. The average H 2 - H 2 interaction energy was then calculated for a given inter- 
molecular distance by taking weighted averages of the interaction energies of these four 
configurations. In five of the configurations (fig. 1(b), (d), (e), (g), and (i)) one of the 
H 2 molecules had a nonequilibrium internuclear distance so that the derivatives of the 
dipole moment with respect to internuclear distance could be calculated from the electric 
dipole moments of all nine configurations . The calculations in this report differ from 
those of references 9 and 10 in the configurations calculated and in mathematical details . 

The author is indebted to Professor V. Magnasco of the Istituto di Chimica In- 
dustriale deU'Universita, Genoa, for providing listings and a deck of various computer 
programs for evaluating molecular integrals. Although these programs were not used 
for this report, they were useful for checking and provided ideas for the programs that 
were used. 


ANALYSIS 

Model and Configurations 

The interaction of two H 2 molecules may be calculated by the orthogonalized 
valence -bond theory of McWeeny (refs. 23 and 24), which utilizes symmetrically ortho- 
gonalized atomic orbitals in the Born-Oppenheimer approximation. In this method 
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orthogonalized orbitals are formed from linear combinations of atomic orbitals. These 
are multiplied by spin eigenfunctions to give orthogonalized spin-orbitals. Formal 
orthogonalized valence -bond structures are formed from linear combinations of anti- 
symmetrized products of the orthogonalized spin-orbitals. The system wave function 
consists of linear combinations of these structures. This involved procedure merely 
amounts to taking the proper linear combination of the octuple products of four atomic 
orbitals and four spin eigenfunctions. This method has the advantage that it allows full 
configuration interaction without serious complexity. It has previously been applied to 
the interaction of two H 2 molecules by Magnasco and Musso (refs. 8 to 10). 

The simplest possible model for a full valence-bond calculation was chosen. The 
two H 2 molecules were assumed to be in the ground electronic state so that for infinite 
separation each of the H 2 molecules could be described by the covalent -ionic valence- 
bond model of Weinbaum (ref. 25). In Weinbaum' s model, two Is atomic orbitals with 
equal orbital exponents (screening constants) were used. To be consistent, in the full 


valence -bond model used for H 2 - H 2 a set of four Is atomic orbitals Xp 


centered on 


the protons was selected. The respective orbital exponents if were assumed to be the 

ir 

same as in the respective isolated H 2 molecules described by the Weinbaum model. If 
both molecules had the same internuclear distances, all orbital exponents were equal. 

If the two molecules had different internuclear distances, two orbital exponents had one 
value, and the other two had a different value. The four normalized Is atomic orbitals 
Xp were 


/? 3\ 1//2 .{r . 

Xp(i) = (-2-J e P P1 p = a, b, c, d (1) 

where the i in parentheses indicates the i th electron, and r pi is the distance between 
proton p and electron i (symbols are given in appendix A). An orbital Xp(i), of course, 
reduces to the wave function for the ground state of an H atom if £ =1.0 bohr~*. 

Four equilibrium configurations (fig. 1(a), (c), (f), and (h)) were chosen so that 
interaction energy could be conveniently averaged over molecular orientations. The 
experimental equilibrium internuclear distance of 1.401446 bohrs (0.0741599 nm) for 
H 2 (refs. 26 and 27) was used rather than Weinbaum’s value of 1. 4166 bohrs 
(0. 074962 nm) (ref. 25) because of the ultimate application of this report to pressure 
induced absorption. 

The other five configurations in figure 1 possessed a nonequilibrium internuclear 
distance so the derivatives of the electric dipole moment with respect to internuclear 
distance could be calculated from the nine configurations. 

In all nine configurations the intermolecular distance R r was varied. It had values 
of 2. 5, 3. 5, 4. 5, and 5. 5 bohrs (1 bohr = 0. 0529167 nm). 
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Energy and Interaction Energy 


To find the electronic energy of the Hg - Hg complex, a Lowdin transformation was 
performed on the orbitals, resulting in orthogonalized orbitals. Formal orthogonalized 
valence -bond structures were formed, and the matrix elements of the Hamiltonian were 
found. An eigenvalue problem was then solved to obtain the electronic energy. The 
nuclear repulsion energy was added to this to get the energy of the complex. The ener- 
gies of the two infinitely separated H2 molecules were subtracted from the energy of the 
complex to yield the interaction energy. These steps are given in the following para- 
graphs. 

Lowdin transformation. - The symmetrical orthogonalization procedure of Lowdin 
(ref. 28) was applied to the Xp- The overlap integral is 

S pq ’ f XpMXqOMvi ® 


Let 


pq 2 ^ 


(3) 


where R is the distance between protons p and q. If £ = £ , S is given by 

PH r 4 r4 

Slater (ref. 3, p. 50). If f / 


s 

pq 4 


e* - e-« 


^_[l + Wnn + Zpq) +e ^[i + i- + -l 

2 \ m 2 / \q -2 -3, 

w pq \ ' \ q q 1 


+ e q-i + A-A 
* q 2 q 3j 


e" w pq 


(4) 


w 


pq 


where 




(5) 
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- 1/2 

It is desired to find the Lowdin transformation matrix ' such that 


X’=XM' 1/2 (6) 

where x is a row vector with elements Xp> and }C_ is a row vector of the four ortho - 
gonalized orbitals x' • First the eigenvalue problem 

USU = d (7) 


is solved, where S is a matrix with elements S . The matrix d is diagonal, with 
— F4 

eigenvalues as diagonal elements. Thus 


,- 1/2 _ 


.- 1/2 

l ll 

0 

0 

0 


0 

a - 1/2 
*22 

0 

0 


0 

0 

d-i/2 

a 33 


0 

0 

0 

d-V2 

a 44 


( 8 ) 


Finally, 




(9) 


Use of equation (6) then accomplishes the Lowdin transformation to obtain X' • 

Formal orthogonalized valence -bond structures . - To produce the proper antisym- 
metry, the x^j are combined in antisymmetrized products of orthogonalized spin- 
orbitals. 

(pqrs) = (4’.)" 1/2 detjxj,(l)a(l), X^(2)a(2), X^(3)/3(3), x^(4)/3(4)] (10) 

where of(l) is the spin eigenfunction of electron 1 with the component of spin angular 
momentum along the axis of quantization equal to B/2, and /3(3) is the spin eigenfunction 
of electron 3 with the component of spin angular momentum along the axis of quantization 
equal to -K/2. The shorthand notation is given on the left side of equation (10), wherein 
a bar indicates /3 spin and no bar indicates a spin. To facilitate computing for some 
purposes, the (pqrs) are arranged in dictionary order in a 6 by 6 matrix (ref. 29). 
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(abab) 

(abac) 

(abad) 

(abbe) 

(abb3) 

(abc3) 

(acab) 

(acae) 

(acad) 

(aebe) 

(acbd) 

(accd) 

(adab) 

(adac) 

(adad) 

(adbc) 

(adbd) 

(aded) 

(bcab) 

(bcac) 

(bead) 

(bebe) 

(bebd) 

(bcc3) 

(bdab) 

(bdac) 

(bdad) 

(bdbc) 

(bdbd) 

(bded) 

(cdab) 

(cdac) 

(edad) 

(cd5c) 

(cdbd) 

(eded) 


For our purposes we will regard the (pqrs) as a 36-element row vector with the elements 
numbered down the first column in equation (11), then down the second column, etc. 

When numbered in this way, set 

cp i = (pqrs) i = 1, 2, . . 36 (12) 

where p, q, r, and s are given in equation (11). 

The formal orthogonalized valence-bond structures j/a are linear combinations of 
the These i/a are eigenfunctions of the total electron spin with total spin angular 
momentum of 0 and component of the total spin angular momentum along the axis of 
quantization of 0. They have been given by Magnasco and Musso (ref. 9) and in the 
present nomenclature are 


^1 = \ ^11 ‘ v 16 ' v 2l ' ^26^ 

^11 = 2 ^ fa 5 + ^25^ 

^2 = j (<p 6 - ^i 6 - <p 2 i + ^ 31 ) 

^12 = 2 ^ fa 12 + ^32^ 

^3 = 2 ^ io * ^ 2 <P 

1^13 = 2 1//2 (-<p 2 - <p 7 ) 

^4 = 2’ 1 / 2 (-*„ - <p 27 ) 

^14 = 2 ^ ("^30 " ^35^ 

^5 = 2_1/2 -VlJ 

^15 = _<;p 8 

= 2 1 / 2 (-<p 2 3 " ^28^ 

^16 = '^15 

^7 ~ ^ (- 4 - <p^g) 

^17 = '^22 

1/2 

^ 8 = 2 (^18 + ^ 33 ) 

^18 = ”^29 

^g =2 ,// (^?2 + 

^19 = ^1 

^ 10 = 2 ^24 ■ < < p 34^ 

^20 = ^ 36 
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where formal singly polar and doubly polar structures are included. All j/a are ortho- 
normal (see ref. 30) except i//j and These may be made orthonormal by the Schmidt 

orthogonalization procedure, which results in replacing i//g in equation (13) by a new 
i//g such that 


_ ^6 + ^ 16 ^21 [ ^26 t ^31 

2 3 1//2 2 • 3 1 / 2 2 • 3 1//2 2 • 3 1 / 2 2 • 3 1 / 2 3 1//2 

Equations (13) (excluding t/zg) and (14) can be represented by 


(14) 


\p = V 


(15) 


where V is a 36 by 20 matrix. The wave function of the system was taken as 

* = ££ (16) 

where C is a column vector determined by minimizing the electronic energy E g . This 
is accomplished by the variation principle for orthonormal functions (ref. 31) which poses 
an eigenvalue problem 


HC=E e C (17) 

where the smallest of the several eigenvalues E g is the ground state electronic energy, 
and 13 is a matrix whose elements are given in the following sections. 

Matrix elements of the Hamiltonian . - The Hamiltonian operator for the electronic 
energy in atomic units with the hartree as the unit of energy is 


A 

H = 





(18) 


o 

Here V. is the Laplacian operator for the coordinates of electron i, and the p sum- 
mation is over protons a, b, c, and d. The elements of the H matrix in equation (17) 
are then 


H 


nk 



K * K dr 


(19) 
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where r includes electron configuration and spin coordinates. From equations (15) 
and (19) 


H = VH V 
— — (p— 


( 20 ) 


where 


(H ) 

v cp 


13 



( 21 ) 


The matrix elements (H ) 

* ij 

a typographical error) if two 


can be found from 
kinds of integrals 


Slater's rules (ref. 30 after correction of 
are known 



where r ^ is the distance from electron 1 to electron 2, dvj is the volume element for 
electron 1, and dv 2 is the volume element for electron 2. From equations (6), (22), 
and (23) 






0=1 


p=l 


10 


( 25 ) 



where 



G 


mnop 



X m (l)x n (l)X 0 (2)x p (2) 


— — dvjdvg 
r 12 


(28) 


The methods of evaluating K^, -'-‘jkp’ anc * ^mnop are &* ven “ appendixes B to D, re- 
spectively. 

Energies . - The energy E is the electronic energy plus the nuclear repulsion 
energy. 


E = E, 


R 


ab 


R. 


ac 


R 


ad 


R 


be 


R 


bd 


R 


cd 


(29) 


The interaction energy E in ^ is the energy E minus the energies E m (r ^) and E 0^) 
of the isolated Hg molecules with internuclear distances r ^ and rg, respectively. 

E int = E - E m< r l> - E n, (r 2> < 3 °) 


The diatomic internuclear distances used in equation (30) were the same as in the Hg - Hg 
complex (see fig. 1). The Hg energies E m (r) were calculated for the covalent-ionic 
valence -bond model (ref. 25). 


Average Interaction Energy 

The average interaction energy E^ of two Hg molecules is, to a first approxi- 
mation, the interaction energy for equilibrium diatomic internuclear distances averaged 
over all possible orientations of the two Hg molecules. Since only four configurations 
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Side view 


End view Configuration 



Figure 2. - Nine equilibrium orientations with equal probability if 
differences in interaction energy are neglected. These orienta- 
tions were used to determine weighting factors for calculating 
average interaction energy (see fig. 4). Configuration in fig- 
ure 1 to which each of the orientations shown is equivalent is 
given at the right. 


with equilibrium diatomic internuclear distances were considered in this report, an 
additional approximation was necessary. Instead of using all possible orientations, the 
nine orientations in figure 2 were used. These are equally probable if differences in 
interaction energy are neglected, and each is equivalent to configurations A, C, F, or 
H (see fig. 1). Hence, 


E. , = — E. , + — E. , + — E. , f + — E. , , (31 

xnt g mt, a g mt, c g int, f g int, h ' 

where the constants are weighting factors. A better value of E. ^ can be found by in- 
cluding Boltzmann factors in the weighting (ref . 32, pp. 921-922), but then E int would 
be complicated by temperature dependence. 
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Electric Dipole Moment 


Calculation of the dipole moment requires an electron population analysis in terms 
of the original nonorthogonal orbitals Xp- A column vector CM may be defined by 

* = (32) 

From equations (15), (16), and (32) 


C’ = VC 


(33) 


The electron population in the orbital density XpX^ irrespective of spin state is 


36 36 / \ 

C, Z ^ 

r=l s=l ' ' 


(34) 


where n 1 ” 011 ® and n are coefficients of spin-orbital products and are called 
ii if rs i* s 

rs ^y Magnasco and Musso (refs. 8 and 29), who give rules for evaluating them 
(ref. 29). 

The electron population in the orbital density XpX^ irrespective of spin state is 
(ref. 8) 



Z K 1/2 ) Pm Z p “ K 1/2 Z 

m = 1 n=1 


(35) 


The normalized 1 -electron population is 



= P S 

pq pq 


(36) 


where the four diagonal elements of the matrix give the atomic populations and the 
off-diagonal elements give the overlap populations. For instance, the population in the 
overlap between orbitals centered on protons a and b is + ^ a - The elements of 
the ^ matrix, of course, add to 4, the number of electrons. 

To find the dipole moment from it is necessary to know locations of the atomic 
charge centers and the overlap charge centers. Cartesian coordinates with the origin 

it. 

half way between the molecules were chosen (see fig. 1(a)). The j Cartesian coordin- 
ate of either type of charge center was designated C -, where p and q correspond to 

P4J 


13 


I 


the two orbitals Xp and involved in an overlap charge or where p and q are equal and 
correspond to the one orbital involved in an atomic charge. The for atomic charges 

are, of course, identical to the coordinates of the protons C .. For overlap charges 

J*J 




+ 




(37) 


where A is the displacement of the overlap charge center from the p-q geometric 
center, measured in the direction from p towards q and expressed in units of R . 
If S p = 5 q , Apq = ° by symmetry. If ? p t C q 



The j th component of the dipole moment in atomic units is then 


(38) 


Mi 


-Z EVrt 

p=l q=l 


(39) 


By symmetry, only the z component may be nonzero. 


Computer Programs 

Calculation of the energy, interaction energy, and electric dipole moment of H 2 - Hg 
was programmed in FORTRAN IV, IBM version 13, for an IBM 7094 digital computer. 
One program calculated all the three- and four -center molecular integrals and had an 
average running time of 8.09 minutes per case. A second program did the remaining 
calculations and had an average running time of 0. 38 minute per case. 
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RESULTS AND DISCUSSION 


In this section results are given for the covalent -ionic valence -bond model of H 2 - 
These results are then used to calculate H2 - H 2 interaction energy and electric dipole 
moment. Comparisons are made with the results of other investigators. 


Covalent-Ionic Valence-Bond Model of H2 

The treatment of the covalent -ionic valence -bond model of H 2 given by Weinbaum 
(ref. 25) was reformulated so values of £ and E m could be calculated for given inter- 
nuclear distance r. The computer program of Dellepiane, Ferro, and Magnasco 
(ref. 33) was used to calculate the necessary integral logarithms. The results are given 
in table I. For an internuclear distance of 1. 4166 bohrs (0.074962 nm), we got 


TABLE I. - ENERGY AND ORBITAL EXPONENT FOR COVALENT- 


IONIC VALENCE -BOND MODEL OF HYDROGEN MOLECULE 


Internuclear distance, 

Energy, 

Orbital exponent, 


r 


m 



bohrs 

m 

hartrees 

j 

bohrs 1 

-i 

m 

1. 300000 

6.8792X10' 11 

-1.14462 

-4. 9899X10' 18 

1.2240 

2. 3131X10 10 

a l. 401446 

7.4160 

-1. 14779 

-5.0037 

1.2002 

2.2681 

b l. 416600 

7.4962 

-1.14790 

-5.0042 

1. 1970 

2.2620 

1.450000 

7.6729 

-1. 14787 

-5.0041 

1. 1896 

2.2481 

1.530000 

8.0963 

-1. 14646 

-4.9979 

1. 1730 

2.2167 


a Gives minimum energy based on spectroscopic data (ref. 26). 
b Gives minimum energy for this model (ref. 25). 


£ = 1. 197 bohrs - *, but Weinbaum got £ = 1. 193 bohrs - *. However, our energy agrees 
exactly with Weinbaum' s. The energy was not very sensitive to the value of £. Values 
of E m and £ from table I were used in the following H 2 - H 2 calculations. 


H 2 - H 2 Energy and Interaction Energy 

Energy and interaction energy of H 2 - H 2 were calculated from equations (29) and 
(30), respectively, and are given in table n for nine nuclear configurations and four 
intermolecular distances. For a given intermolecular distance, the linear configurations 
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TABLE H. - ENERGY AND ELECTRIC DIPOLE MOMENT FOR TWO HYDROGEN MOLECULES 


Ihtermolecular 

distance, 

R r 

Con- 

figura- 

tion 3 

Energy, 

E 

Interaction 

energy, 

R int 

z -component of 
electric dipole 
moment* 5 

bohrs 

m 


hartree 

J 

hartree 

J 

at. units 

C-m 

2.5 

1.3229X10 -10 

A 

-2.04492 

-8. 9147X10" 18 

0.25066 

1.0927X10" 18 

0 

0 



B 

-2.03489 

-8.8710 

. 26078 

1. 1369 

.034890 

2.958X10" 31 



C 

-2.15035 

-9. 3743 

.14523 

.6331 

-.347786 

-2. 949xl0~ 30 



D 

-2. 15018 

-9. 3736 

. 14549 

.6343 

-.343797 

-2.915X10" 30 



E 

-2. 14658 

-9. 3579 

. 14909 

.6499 

-. 384906 

-3.263X10' 30 



F 

-2.18153 

-9. 5102 

.11406 

.4972 

0 

0 



G 

-2. 18060 

-9.5062 

.11506 

. 5016 

.004599 

3.899X10' 32 



H 

-2. 18434 

-9.5225 

.11124 

. 4849 

0 

0 



I 

-2. 19410 

-9.5650 

.10157 

.4428 

.003720 

3.154X10" 32 

3.5 

1.8521X10' 10 

A 

-2.25667 

-9.8378X10* 18 

0.03891 

0. 1696X10' 18 

0 

0 



B 

-2.25540 

-9.8323 

. 04026 

. 1755 

.010739 

9. 104X10" 32 



C 

-2.26824 

-9.8882 

.02734 

. 1192 

-.097871 

-8.297X10 -31 



D 

-2.26800 

-9.8872 

.02766 

. 1206 

-.097330 

-8.251X10' 31 



E 

-2.26735 

-9.8844 

.02832 

. 1235 

-.109437 

-9.278X10' 31 



F 

-2.27440 

-9.9151 

.02118 

.0923 

0 

0 



G 

-2.27412 

-9.9139 

.02155 

.0939 

.001596 

1.353X10" 32 



H 

-2.27456 

-9.9158 

.02102 

.0916 

0 

0 



I 

-2.27706 

-9.9267 

.01861 

.0811 

.001287 

1.091X10* 32 

4.5 

2.3813X10' 10 

A 

-2.28895 

-9.9785X10' 18 

0.00663 

0.0289X10" 18 

0 

0 



B 

-2.28874 

-9.9776 

. 00692 

.0302 

. 002468 

2.092X10" 32 



C 

-2.29084 

-9.9868 

. 00475 

.0207 

-.024428 

-2.071X10" 31 



D 

-2.29082 

-9.9867 

. 00485 

.0211 

-.024457 

-2.073X10' 31 



E 

-2.29070 

-9.9861 

.00497 

.0217 

-.027478 

-2. 330X10' 31 



F 

-2.29194 

-9.9916 

.00364 

.0159 

0 

0 



G 

-2.29193 

-9.9915 

.00374 

.0163 

.000457 

3.87X10' 33 



H 

-2.29193 

-9.9915 

.00365 

.0159 

0 

0 



I 

-2.29257 

-9.9943 

.00310 

.0135 

.000357 

3.03X10" 33 

5.5 

2.9104X10" 10 

A 

-2.29453 

-10.0028X10" 18 

0.00105 

0. 0046X10” 18 

0 

0 



B 

-2.29456 

-10.0030 

.00111 

.0048 

. 000498 

4.22X10' 33 



C 

-2.29484 

-10.0042 

.00074 

.0032 

-.005873 

-4.979X10" 32 



D 

-2.29490 

-10.0045 

.00076 

.0033 

-.005938 

-5.034X10" 32 



E 

-2.29488 

-10.0044 

.00078 

.0034 

-.006616 

-5.609X10' 32 



F 

-2.29500 

-10.0049 

.00058 

.0025 

0 

0 



G 

-2.29506 

-10.0052 

.00060 

.0026 

.000111 

9.4X10" 34 



H 

-2.29499 

-10.0048 

.00059 

.0026 

0 

0 



I 

-2.29519 

-10.0057 

. 00048 

.0021 

. 000084 

7. 1X10' 34 


a See fig. 1. 

^See fig. 1(a) for coordinates. 


Other components of the electric dipole moment are zero by symmetry. 
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Interaction energy, Ej n p J Interaction energy, Ej n p J 


I 


Source of theory 

O Taylor (ref. 7) 

□ Griffing and Vanderslice (ref. 4) 

This report <eq. (30)) 

Schwartz and Schaad (ref. 14) 

Ruffa and Griffing (ref. 6) 

Mason and Hirschfelder (ref. 17) 


^ \ 


v \ 

X\ 


(a) Linear configuration A (see fig. 1(a)). 


This report (eq. (30)) 

Mason and Hirschfelder (ref. 17) 


i 1 i I — 

3. 5 4.0 4.5 5.0 

Intermolecular distance, R r , bohr 


1.5 20 2.5 3.0x10 

Intermolecular distance, R r , m 

(b) Planar perpendicular configuration C (see fig. 1(c)). 

Figure 3. - Interaction energy for four configurations of two hydrogen molecules, each 
configuration having equilibrium internuclear distances in both molecules. 



Interaction energy, E int , J Interaction energy, E: 








had the highest interaction energies, followed by the planar perpendicular configurations 
The quadrilateral and nonplanar perpendicular configurations had the lowest interaction 
energies. It is not safe to draw any conclusions about the change in interaction energy 
between equilibrium and nonequilibrium configurations because the experimental equili- 
brium internuclear distance was used in this report rather than Weinbaum's value 
(ref. 25 and table I). Values of interaction energy in table II are plotted in figure 3 for 
the four equilibrium configurations. The interaction energies varied approximately 
exponentially with inter molecular distance. 


H2 - H2 Average Interaction Energy 

The average interaction energy of Hg - Hg was calculated from equation (31) for four 
intermolecular distances. Results are given in table HI and are plotted in figure 4. 


TABLE m. - AVERAGE INTERACTION 
ENERGY FOR TWO HYDROGEN MOLE- 
CULES WITH EQUILIBRIUM INTER- 
NUCLEAR DISTANCES 11 


Intermolecular 

Average 

distance, 

interaction energy, 


R r 

i 

^int 

bohrs 

m 

hartrees 

j 

2.5 

1. 3229X10* 10 

0. 14246 

6.2ioxio' 19 

3. 5 

1. 8521 

.02585 

1. 127 

4. 5 

2.3813 

.00447 

. 195 

5. 5 

2.9104 

.00071 

.031 


Equilibrium internuclear distances of both 
molecules were 1.401446 bohrs 
(7. 4160X10' 11 m) (ref. 26). 


H 2 - H 2 Electric Dipole Moment and Its Derivatives 

The electric dipole moment of Hg - H2 was calculated from equation (39) for the 
nine configurations in figure 1 and four intermolecular distances. The x and y com- 
ponents were all essentially zero, as required by symmetry. The z component is given 
in table n. The planar perpendicular configurations had the largest magnitude of dipole 
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Intermolecular distance, R r , m 


Figure 4. - Average interaction energy for two hydrogen molecules. The curve labelled 
"This report" was obtained by averaging interaction energies for the nine equally 
probable orientations given in figure 2. 


moment because of their high degree of asymmetry. Equilibrium configurations A, F, 
and H had zero dipole moment by symmetry. The z component of dipole moment of 
equilibrium configuration C is plotted in figure 5. It varies approximately exponentially 
with intermolecular distance. The derivatives of the dipole moment with respect to 
internuclear distances were calculated from the dipole moments of the nine configurations 
(see table II) and are shown in figure 6. 


Comparison of H2 - H2 Results with Other investigators 

Interaction energy . - Comparisons of from equation (30) with ab initio results 

of other investigators are given in figure 3 for four H2 - Hg configurations with equil- 
ibrium internuclear distances in both molecules. 
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This report (eq. (39)) 

Large R r approximation using Q s 

from ref. 34 and a from ref. 35 



Intermolecular distance, R r , bohrs 


1.5 ZO 2J5 i Ox 10' 10 

Intermolecular distance, R r , m 

Figure 5. - z component of electric dipole moment for planar perpendicular configuration C 
(see fig. 1(c)) of two hydrogen molecules. Coordinates are shown in figure 1(a). 


Figure 3(a) is for the linear configuration. The interaction energies of Griffing and 
Vanderslice (ref. 4) and Ruffa and Griffing (ref. 6) are high, possibly because of in- 
accurate molecular integrals. The energy of Taylor (ref. 7) is also high, possibly be- 
cause he set all the orbital exponents equal to one. The energies of Schwartz and Schaad 
(ref. 14) are lower than our energies. The variation principle states that E is an upper 
limit, but gives no information on E int - Consequently, this principle cannot be used to 
show that Schwartz and Schaad's is best. However, since they used 24 orbitals 

compared with four in this report, their E int is probably best for the range of R r they 
covered. The theory of Mason and Hirschfelder (ref. 17) is believed best at large R f 
but is not accurate at small R r because it is principally a perturbation calculation. Its 
superiority at large R r is due to its accurate semiempirical treatment of dispersion 
energy. 

Figures 3(b) to (d) are for the planar perpendicular, quadrilateral, and nonplanar 
perpendicular configurations, respectively. The E^ of Mason and Hirschfelder 
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Intermolecular distance, R r , bohrs 


1.5 ZO a 5 3. 0x10 ' 10 

Intermolecular distance, R r , m 

Figure 6 . - Derivatives of electric dipole moment with respect to internuclear distance for 
various configurations of two hydrogen molecules. Internuclear distances rj and r 3 
are the distances between the nuclei of molecules 1 and 2, respectively (see fig. 1 ). 

(ref. 17) are believed best at large R p . For the quadrilateral configuration there is no 
discernible difference between Magnasco and Musso's (ref. 9) and our E int . 

Average interaction energy . - Comparison of Ej t from equation (31) with results 
of other investigators are given in figure 4 for Hg - Hg with equilibrium internuclear 
distances in both molecules. The Lennard- Jones 6-12 potential (ref. 32) was fitted to 
data from moderate temperature equation-of -state experiments and hence is not accurate 
at small R r . The Morse potential (ref. 22) was fitted to moderate temperature equation- 
of-state and viscosity data and also to viscosity data up to 2340 K (ref. 21); therefore, it 
should be valid down to R r = 4. 0 bohrs (2. 1x10"^ m). However, it is obviously wrong 
for small enough R p because it does not approach infinity as R f approaches zero. 

Thus its agreement with equation (31) for 2. 5<R r < 4 bohrs (0. 13== R r ^ 0.21 nm) is 
surprising. Two semiempirical theories (refs. 17 and 18) are included in figure 4 for 
reference only. I believe the excessive values of E- nt from equation (31) for large R r 
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are due to insufficient allowance for electron correlation in the wave function. This 
causes the magnitude of the dispersion energy (which is a negative quantity) to be low. 

Electric dipole moment. - Comparisons of a from equation (39) with results of 
other investigators are given in figure 5 for the planar perpendicular configuration. The 
large R r approximation was obtained by considering the polarization of each molecule 
by the electric quadrupole of the other. The scalar quadrupole moment Q„ and polar- 

b 

izability a of H 0 were obtained from Kolos and Wolniewicz (refs. 34 and 35) and inspire 
“ 4 

confidence. The resulting p z varies as 1/R r as shown and is correct without much 
doubt for large enough R r . A curve is also given in figure 5 for the large R r approxi- 
mation plus a semiempirical overlap contribution (ref. 20). The overlap contribution is 
not large, so the trend is essentially the same as for the large R r approximation alone. 
Equation (39) does not follow either dashed curve for R r >4. 5 bohrs (0. 24 nm) and, in 
fact, is even curved the wrong way. I believe this difficulty is due to the model used: 
at large R r it essentially describes, two covalent-ionic valence-bond Hg molecules, 
which have no polarizability perpendicular to their axes and probably have incorrect 
quadrupole moments. (The quadrupole moment is very sensitive to the electronic charge 
distribution.) Consequently, the orthogonalized valence-bond theory used in this report 

does not give reliable a for R >4.5 bohrs (0.24 nm) despite the inclusion of singly 

z r 

and doubly polar structures. 


CONCLUDING REMARKS 

Full valence -bond ab initio calculations of Hg - Hg interaction energy and electric 
dipole moment were carried out with accurate values for all molecular integrals. Lin- 
ear, planar perpendicular, quadrilateral, and nonplanar perpendicular configurations 
were included with inter molecular distances from 2.5 to 5.5 bohr (0. 13 to 0. 29 nm). A 
weighted average of the interaction energies of the four equilibrium configurations was 
then taken. 

The interaction energies obtained appear to be valid for intermolecular distances 
less than 4. 5 bohr (0. 24 nm) but are too large at appreciably greater intermolecular 
distances because of insufficient allowance for electron correlation in the wave function. 

Fisher (ref. 22) has fitted a Morse potential to viscosity and equation-of -state data, 
and this potential should be valid down to intermolecular distances of 4. 0 bohrs (0.21 
nm), based on the temperatures at which the data were taken. The calculations in this 
report show that this Morse potential is also a good approximation for intermolecular 
distances between 2. 5 and 4. 0 bohrs (0. 13 and 0. 21 nm). Consequently, this Morse po- 
tential is useful for practical applications for all R r greater than 2. 5 bohrs (0. 13 nm). 

The electric dipole moments obtained by the full valence -bond calculation appear to 
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be approximately correct for inter molecular distances up to 4. 5 bohrs (0. 24 nm). For 
appreciably larger distances they are too small due to limitations of the model. The 
derivatives of the dipole moment with respect to the internuclear distances of the 
molecules were also obtained. These are necessary to calculate pressure-induced vi- 
brational absorption coefficients . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 23, 1969, 

122-28. 
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APPENDIX A 


SYMBOLS 


A 

definite integral involving ex- 



ponential 


(abed) 

shorthand notation for anti- 



symmetrized product of ortho- 
gonalized spin-orbitals 

G mnop 

C 

column vector of coefficients of 

H 

ci 

column vector of coefficients of 

A 

H 

c p] 

Vi 

11. 

j Cartesian coordinate of 

H k 

proton p 


C PPj 

j ul Cartesian coordinate of 


center of atomic charge p 

H nk 

c pqj 

iU 

j m Cartesian coordinate of 

5* 

center of overlap charge pq 

c’ 

r 

element of C 1 

%).. 
* ij 

d 

matrix of overlap eigenvalues 

R 

d-l/2 

inverse square root of d 

J 

d-V 2 

-1/2 

diagonal element of d 7 

K 

E 

electronic and nuclear repulsion 
energy of complex 

K k 

E e 

electronic energy 


E int 

interaction energy 

K k* 

E int 

average interaction energy 

L jkp 

E m 

electronic and nuclear repulsion 
energy of diatomic molecule 

L i 


one -electron integral of ortho - 
gonalized orbitals 
two -electron integral of ortho - 

M 

G ijkZ 

gonalized orbitals 



definite integral involving ex- 
ponential and Legendre func- 
tion of first kind 

two-electron integral of Is 
orbitals 

Hamiltonian matrix with basis 
functions i/a 

Hamiltonian operator for elec- 
trons 

definite integral involving ex- 
ponential and Legendre func- 
tions of first and second kinds 

element of H 

Hamiltonian matrix with basis 
functions cp ^ 

element of H „ 

_£ 

Planck constant divided by 2n 

center of orbital Xj 

center of orbital 

definite integral involving ex- 
ponential and Legendre func- 
tion of first kind 

kinetic energy integral 

nuclear attraction integral 

definite integral involving ex- 
ponential and Legendre 
function of second kind 

center of orbital x^, 
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M- 1 / 2 

X 

Lowdin transformation matrix 


element of 

V * / ki 

X 

N 

center of orbital X,^ 

„mana 

n rs 

coefficients of spin-orbital 
products with two a spins 

m/3n/3 

rs 

coefficients of spin-orbital 
products with two /3 spins 

P 

center of proton p 


matrix of 

pq 

P k 

Legendre function of first 
kind 

P mn 

electron population in the 
orbital density X^X^ 
irrespective of spin state 

P pq 

electron population in the 
orbital density XpX^ ir- 
respective of spin state 

^pq 

normalized one -electron 
population 


constants in Legendre func- 
tion of first kind 

Q k 

Legendre function of second 
kind 

Q s 

scalar quadrupole moment 

q 

(£p - £ q ) R pq/2 (subscripts 
vary) 

A 

q 

0 

a 

a 

/^S 

O 

1 

S 

^J3 

q 

- V R mp /2 

R pq 

distance between points p 
and q 

R r 

intermolecular distance 
(see fig. 1) 


r ’ r i 

internuclear distance of diatomic 
molecule i 

r 

+ V R mo 

r 

< 3 ?m + VW 2 

r ij 

distance between electron or 
proton i and electron j 

S 

matrix of S 

— 

pq 

V 

overlap integral for Xp and Xq 

T 

definite integral involving ex- 
ponential 

U 

unitary matrix 

u 

tranpose of U 

U mn< 2 > 

electrostatic potential at electron 
2 due to the bicentric charge 
distribution X m (l)x n (l) 

u 

variable of integration 

V 

matrix relating xp and cp 

V 

transpose of V 

v i 

volume in configuration space for 
electron i 

W ’V 

V W 2 

X k’ Z k 

x and z coordinates of point K 
in units of R^-/2 

x,y,z 

Cartesian coordinates 

a 

polarizability 

«(i) 

spin eigenfunction of electron i 
with component of spin angular 
momentum along the axis of 
quantization equal to ft/2 

/3(i) 

spin eigenfunction of electron i 
with component of spin angular 
momentum along the axis of 
quantization equal to -fi/2 
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I 


A pq 

displacement of overlap charge 
center for p-q geometric 
center, measured in direction 
from p towards q and ex- 
pressed in units of R 

Jr4 

? 

orbital exponent 


orbital exponent of Xp 

9, i u,v 

spheroidal coordinates 

\ 

-i/m 


l± of point K 

Moo 

finite upper limit for i± inte- 
gration 

M(0 

smaller of /ij and jUg 

m(>) 

larger of j Uj and /u 2 

T 

volume in configuration and spin 


space of four electrons 


£ 

row vector with elements cp^ 


antisymmetrized product of 
orthogonalized spin-orbitals 

X 

row vector with elements Xp 

Xl 

row vector with elements x^ 

Xpft) 

Is orbital of electron i centered 
on proton p 

X^(i) 

orthogonalized orbital of electron 
i 

* 

wave function of system 


row vector with elements \}/. 


formal orthogonalized valence- 


bond structure 
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APPENDIX B 


EVALUATION OF KINETIC ENERGY INTEGRALS 


Herein methods are given for evaluating kinetic energy integrals involving Is atomic 
orbitals (eq. (1)) with equal or unequal orbital exponents. These integrals all have the 
form 


^ ■/ v»(-i 


Equal Orbital Exponents 

Analytic expressions for for the two cases k = l and k ^ l have been given 
in reference 3 (p. 50) (but a typographical error must be corrected in the k ^ l case). 


Unequal Orbital Exponents 

If k ^ l and the orbital exponents are unequal, the integral in equation (Bl) may be 
found analytically by the methods in reference 3. 



where w^ is given by equation (3) and where 


28 



(B3) 


5 - fa - 

2 

All K kZ were evaluated to within ±1 in the sixth significant digit. 
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APPENDIX C 


EVALUATION OF NUCLEAR ATTRACTION INTEGRALS 

Methods for evaluating nuclear attraction integrals involving Is atomic orbitals 
(eq. (1)) with equal or unequal orbital exponents are given herein. These integrals all 
have the form 


L Jkp- f 

where r p ^ is the distance from proton p to electron 1. 


(Cl) 


Equal Orbital Exponents 

If all three subscripts of Lj kp are equal or if only two of them are equal, analytic 
expressions for L^ kp are given in reference 3 (p. 50). If all subscripts are different, 
numerical integration is required and is described in the following paragraphs. 

Numerical integration of three -center integrals is best done in spheroidal coor- 
dinates (ref. 3). First, however, a set of Cartesian coordinates is selected with the 
origin half way between P and J and with the z axis in the direction from P to J. 


Z 
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The x axis is chosen so K is in the x-z plane with x^ zero or positive. The three 
spheroidal coordinates 9, n, and v can now be defined. The azimuthal angle measured 
around the z axis as shown is called 9 . The other two coordinates are 




f pl + f jl 



v 




(C2) 

(C3) 


where r^j is the distance from P to dvp and r^ is the distance from J to dvj. 
Surfaces of constant ju are ellipsoids of revolution about the z axis with foci at P and 
J. Surfaces of constant v are hyperboloids of revolution about the z axis with foci at 
P and J except for v = 0, which is the x-y plane. 

In these coordinates, equation (Cl) becomes 


L. k =-^l 

i k e 2* 


/»°° /»1 

/ f 


-t 

(m - v)e ^5 2 


jf exp [" -f (*£ 


+ z£ - 1 + /I 2 + v 2 - 2iivZ 


k- 2X k 


1 u 2 - 1 -^1 - v 2 cos 9 


■) 


1/2 


dQdvdy. (C4) 


where and are and z^, respectively, in units of Rpj/2- The 9 integration 
was done by a single Gaussian quadrature . The v integration was done by using two 
Gaussian quadratures: one from -1 to ^ and one from to 1. The ju integration 
was also split into two parts. Integration from 1 to ^ was done by Gaussian quadra- 
ture. From ju^ to an upper limit the integration was done by transforming the 
integral so it was with respect to a variable X = -1/ (i and applying Gaussian quadrature. 
(The author is indebted to Professor V. Magnasco, Universita di Genova, Genoa, for 
this idea. ) In each case an appropriate finite i± x had to be calculated to replace 00 as 
the upper limit in the fi integration. 

Use of the coordinate system shown was superior to use of the coordinate system in 
reference 36 because it gave better accuracy for the same number of integration points 
for cases where J, K, and P were in line or almost in line. This was apparently be- 
cause the integrand in equation (C4) is everywhere finite, which is not the case for the 
equation in different coordinates in reference 36. 


31 



Unequal Orbital Exponents 


If j = k ^ p, the only two orbital exponents that occur in the integral are identical. 
However, if p = j^k or p = k ^ j, the unequal orbital exponents result in new express- 
ions for Ljkp- For the case p = k ^ j 


X 1 Xj2 03 

L ikp = -2 R ikyfj?k 


e q- e -ql + w jk r -qV 1 t 1 


w. 


jk 


U q 2i 


-e^fl -- 1 

^ q 2 1 


■ W jk 


W jk 


(C5) 


where the methods of reference 3 were used and where 


q = 


ft .i - g k> R ik 


(C6) 


For the case where none of the subscripts of L.. are equal, equation (C4) is applicable. 
All were evaluated to within ±3 in the sixth significant digit. 
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APPENDIX D 


EVALUATION OF TWO-ELECTRON INTEGRALS 


Methods are given for evaluating two-electron integrals involving Is atomic orbitals 
(eq. (1)) with equal or unequal orbital exponents. These integrals all have the form 


mnop 


=jj ^nWXnWXo^V 2 ) ~ dv l dv 2 


(Dl) 


where r ^ is the distance between the two electrons. 


Equal Orbital Exponents 

One -center integral . - An expression for when all the subscripts are equal 

is given in reference 3 (p. 65). 

Two- center Coulomb integral . - If m = n and o = p, an expression for G mn0 p is 
given in reference 3 (p. 50). 

Two -center exch ange integral . - If m = o and n = p, an expression for G mnQ p is 
given in reference 3 (p. 50), but it requires correction of a typographical error. The 
expression involves the integral logarithm, which was evaluated using the program of 
Dellepiane, Ferro, and Magnasco (ref. 33). 

Two-center hybr id integral . - An expression for G mn0 p if three of the subscripts 
are identical is given in reference 3 (p. 65). 

Three - center C oulo m b integral . - If m = n / o / p or m / n / o = p, the integral 
must be integrated numerically. The method of Magnasco and Dellepiane (ref. 36) was 
used except that on the outer p integration range the p integral was transformed so 
that it was with respect to a variable X = -1/p and Gaussian quadrature was applied. 

Three-center exchange integral . - If one of the first two subscripts equalled one of 
the last two subscripts of G mn0 p, the integral was integrated numerically by the method 
of Magnasco and Dellepiane (ref. 37) except that on the outer p integration range the 
method of the preceding paragraph was used. 

Fo ur-center integral . - For cases where all the subscripts of G mnQ p were differ- 
ent, an evaluation method has been given by Magnasco and Dellepiane (ref. 37). This 
method was used with three modifications to improve the accuracy: (1) on the outer p 
integration range the p integral was transformed so it was with respect to a variable 
X = -1/p, and Gaussian quadrature was applied; (2) in some cases three p integration 
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ranges were used instead of two; and (3) in some cases three u integration ranges were 
used instead of two. 


Unequal Orbital Exponents 


Two-center Coulomb integral . - Here m = n and o = p. Using the methods in 
reference 3 gives 
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(D2) 


where 


f * «m + Co> R mo 


(D3) 


and 


4 = <C m - 5 0 )E m0 


(D4) 


Two -center exchange integral . - Here m = o and n = p. This can be expressed as 
an infinite series. Spheroidal coordinates were used for both electrons. The quantity 
l/rj 2 was expanded in the Newmann expansion (ref. 3, p. 266) with the result 


G mnop4 R ^ ? "2 (a t l)k(® 2 H k (2,2,wJ 

8 k=0 


2G k (q)G k (q)H k (0 ’ 2 ’ w mn ) + G k (q)2H k (0 ’°’ w mn ) ] < D5 ) 
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where w 


is given by equation (3) and 



Here the P k and Q k are Legendre functions of the first and second kinds, respectively, 
p(<) is the smaller of pj and j u 2 , and p(>) is the larger of pj and / u 2 - The G k 
integrals were evaluated by Rosen (ref. 38). The H k integrals were evaluated by James 
and Coolidge (ref. 39) and Rosen (ref. 40). 

Two-center hybrid integral. - Here, three of the subscripts are identical. This 
case may be integrated by methods in reference 3. If the first three subscripts are 
identical, 



where w mp is given by equation (3) and 



(DIO) 
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(D12) 


Three -center Coulomb integral . - Here, m = n^o^p or m^n^o = p in 
G mn0 p- The method of evaluation was the same as for equal orbital exponents. 

Three-center exchange and four-center integrals . - The methods of evaluating these 
integrals were the same as for equal orbital exponents with the following exception. The 
methods of Magnasco and Dellepiane (ref. 37) require evaluation of 


D nm< 2 > = f X^DVD-i-dV! (D13) 

J r 12 

as a step in their procedure, where U mn (2) is the electrostatic potential at electron 2 
due to the bicentric charge distribution (l)x n (l)- They give an expression for 

U mn (2) u ^m = ? n’ but not ? m ^ ? n A solution for C m t £ n wiU now be given. 

An infinite series for U mn (2) with £ m ^ £ n may be derived by using spheroidal 
coordinates for both electrons 1 and 2. These coordinates were the same as in appen- 
dix C except that M replaced P and N replaced J. The quantity l/r^g w as ex- 
pressed as a Neumann expansion. The result is 


CO 
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and 
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-CC - W|Lt 

, k (w,M 2 ) = / e 

•'Mo 


(D16) 


The integrals may easily be evaluated explicitly. To do this, the Legendre 
function of the first kind may be expressed as 


p k ( ^ = E 4^ 

3=0 

where the p^ are constants. Substituting equation (D17) into (D15) gives 


(D17) 


k 

K k (w ,/j) = Yj I 4a(w,juJ + i) 
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where 


A(w }J u, 0) = i (-e" w/i + e” w ) 
w 


(D 18 ) 


(D19) 


and the recursion relation is 


A(w,/u,h) = — £-/i k e“ w ^ + e" w + hA(w, ju, h - lj] 
w 


(D20) 


The integrals are not so easy to evaluate to six significant digits . The Legendre 
function of the second kind may be expanded in a hyper geometric series 


Q k (M) = 
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M k+1 (2k + 1)1 » 
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which is abbreviated as 



(D22) 


Substituting equation (D22) into (D16) gives 




ki 

(2k + 1)!’. 




(D23) 


where 


T(w,ju,h) s / e -wu u h du (D24) 

and h can be a positive or negative integer or zero. Equation (D24) was integrated 
analytically to get T in terms of exponentials and integral logarithms. If 1. 2, 
equation (D23) gives six significant digits without too many terms. For 1. 2, it is 
better to use 



e" wu Q k (u)u 1 du + l£(w, 1.2) 


(D25) 


where the first and second terms on the right were evaluated by Gaussian quadrature and 
equation (D23), respectively. 

General . - All G mnop were evaluated to within ±4 in the sixth significant digit. 
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